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ABSTRACT 

During a supernova explosion, a radiation-dominated shock (RDS) travels through 
its progenitor. A collisionless shock (CS) is usually assumed to replace it during 
shock breakout (SB). We demonstrate here that for some realistic progenitors en¬ 
shrouded in optically thick winds, such as possibly SN 2008D, a CS forms deep inside 
the wind, soon after the RDS leaves the core, and therefore significantly before SB. 
The RDS does not survive the transition from the core to the thick wind when the 
wind close to the core is not sufficiently dense to compensate for the r -2 dilution 
of photons due to shock curvature. This typically happens when the shock velocity 

is ~ °- lc ( i 0 km/s )( 5-10-4M 0 /yr )~ 1 ( i O iAm )» where u w, M and r* are respectively the 
wind velocity, mass-loss rate and radius of the progenitor star. The radiative CS re¬ 
sults in a hard spectrum of the photon flash at breakout, which would produce an 
X-ray flash. Cosmic ray acceleration would start before SB, for such progenitors. A 
fraction of secondary TeV neutrinos can reach the observer up to more than ten hours 
before the first photons from breakout, providing information on the invisible layers 
of the progenitor. 

Key words: acceleration of particles - plasmas - shock waves - cosmic rays - su¬ 
pernovae: general. 


1 INTRODUCTION 

Type Ib/c and If supernovae (SNe) are generated by core 
collapse in massive stars. When the central engine forms, a 
shock wave is launched through the hydrostatic core of the 
progenitor. The shock is radiation-dominated (or radiation- 
mediated), i.e. the radiat ion pressure in the down stream 
exceeds the fluid pressure rtZel’dovich fe Raizeij|l966h . Once 
the radiation-dominated shock (RDS) reaches the optically 
thin outer layers of the stellar core or of its wind (if opti¬ 
cally thick), photons cannot stay confined in the immediate 
downstream and escape ahead of the shock. This flash of 


Falk! Il978l: Klein & Chevalier 

19781: Chevalier & Klein 

19791: Ensman & Burrows 

1991: Matzner & McKee 

19991; Blinnikov et al. 200d: Calzavara & Matzner 20041: 

Waxman et al.j2007: Katz et al. 

2014 20121: Piro et al. 2010i; 

Nakar & Sari 2010: SaDir et al. 

l201ll. 20131). Up until now 


a few of them hav e been observe d (tCjmrBana^et_alJ |2006|; 


Gezarigtjd, 2Q0fl; ModjazeLah 2009|; Schawinski et al 


20081 : ISoderberg et al.ll2008l : lOfek et al.l 20101b and some X - 
ray fla s hes (XRFs) may b e related to SB ( e.g. KufkamU^il 
d 1998lf ; I Tan et ahl (1200 il l ; iMazzali et al.l d2008ll : Katz et al 


d201lll l. See lOfek et al.l (l2013al lbh and iMurase et all d2014l f 
for radiative signatures at and following breakout. 

At SB, the RDS disap pears and a co llision¬ 
less shock (CS) l ate r for ms ( Ch evalje^jfeJKjeinl 


Ensman fe Burrowsi 19921: Iwi 


fe Loebl 


197S: 


2001 


Chevalier fe Franssonl l2008l l. The Larmor radius 7 ~l of 


suprathermal particles is smaller than the width of the 
RDS, which is ~ \ c/3u s for a sh ock velocity u s and photon 
mean free path A dWeavei Il976lf. On the other hand, tl is 
larger than the CS width ( Belli 1978al fbh and diffusive shock 
acceleration becomes possible. A thorough understanding 
of the CS formation time is then crucial to study the 
onset of CR accelerat ion, when very high en e rgies might 
be reached: > TeV dWaxman fe Loebl j200ll: iKatz et aid 
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1201 lTl . PeV dTatischefll 20091 : Bell et ail] 120131') . and mayb e 
ultra-high energies for transrelativistic SNe ( Budnikll2008l f. 

Post-main-sequence mass-loss of massive stars is suffi¬ 
ciently high for some SN progenitors, such as some Wolf- 
Rayet (WR) stars, blue and red supergiants ( RSG), to 
end u p surrounded with optically thick winds dCrowtheid 
120071 : lLangeil [2012). Also, remar kable outbursts ca n oc- 
cur before the explosi on, see e.g. lOfek et ~ahl d2013bfl and 
ISvirski fe N akail d20l4l . For optically thick winds, the hy¬ 
drostatic surface is not observable, which complicates our 
understanding of the late stages of massive star evolu- 
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tion dGroh et al.l l2013h . iKatz et al.l d201ll ') demonstrated 
that a CS must appear d uring or on t h e time scale of SB (se e 
also ICh evalier fc Irwinl (12011. l2012f). iMurase et al.l l|201 il l. 
ISvirski et all (l2012h and lKashivama et al.l (l2() IOl ) ). when the 
RDS reaches the optically thin layers of the win d, at an op¬ 
tical depth r ~ c/u s = (&T 1 from the surface dKatz et al.1 

l20ldl . 

We demonstrate in this paper that the formation of a 
collisionless shock occurs significantly before SB for some 
progenitors enshrouded in optically thick winds: For some 
realistic density profiles and shock velocities, the RDS stalls 
in the optically thick layers of the wind. A CS forms within 
the broader radiation-dominated transition soon after the 
RDS leaves the hydrostatic core of the progenitor, at r > 
/3s 1 - This makes SB from some optically thick winds funda¬ 
mentally different from thin winds. From a theoretical per¬ 
spective, it redefines the onset of CR acceleration with re¬ 
spect to SB, since it can start in such cases significantly 
before SB. This provides a new method to constrain ob- 
servationally otherwise inaccessible parameters such as the 
radius of the invisible hydrostatic core, see below. Also, the 
spectrum and energy emitted from the beginning of SB are 
affected by the earlier presence of a radiative collisionless 
shock and secondaries of 7-rays from CRs injected deep 
within the wind. 

Supernova SN 2008D/XRF 080109 may have been an 
event in which a CS is formed befor e SB, assuming a progen - 
itor with the parameters derived in lSvirski &; Nakarl (l2014|j . 

The paper is organized as follows. In Section [2j we re¬ 
call a few results for progenitors surrounded with optically 
thin winds. In Section]]!] we discuss conditions for the forma¬ 
tion of a collisionless shock before breakout from optically 
thick winds, and confirm these with numerical simulations 
in Section [T] We investigate in Section [5] particle accelera¬ 
tion in thick winds, and discuss in Section [6] observational 
consequences of our findings. 


2 OPTICALLY THIN WINDS 

In the following, we assume a non-relativistic shock in spher¬ 
ical symmetry with radius r, where r = 0 corresponds to the 
centre of the progenitor. The hydrostatic core and the wind 
are assumed to be fully ionized hydrogen. Assuming more re¬ 
alistic compositions would not change our findings. For the 
temperatures we consider, between ~eV and ~ m e c 2 , the 
opacity k is dominated by Thomson scattering: k = at/m p , 
where cr t is the Thomson cross section. 

In this Section, we first consider a progenitor sur¬ 
rounded with an optically thin wind. SB then starts in the 
outer layers of the core at r « c/3m s dKlein fc Chevalier! 
1 19781 : [Chevalier fc Kleinll l970t) . For such progenitors, CS are 
not expected to form before SB. 

I11 Lagrangian coordinates, the acceleration of a shell 
of wind is Du/Dt = nT^d/c — (1/p) dp/dr, where _F ra d is 
the photon flux, u the shell velocity and p the fluid pres¬ 
sure. The maximum velocity that can be reached by a shell, 
initially at n, due to the flash of photons from breakout 
at thi is Umax,7 = k f t °° -Frad dt/c < k Cdt/Ancrf oc rf 2 
(see also IKatz et al.1 (folll B where C denotes the SN lu¬ 
minosity. After beginning of SB at r^r, the formation of a 
collisionless shock is not immediate, see [Klein fe Chevalier! 


(1 1978l l and IChevalier fe Franssonl J2008l l. The r~ 2 dilution 
of breakout photons ensures that a shell <Si initially at 
ri sS rv will catch up supersonically a shell S 2 initially at 
r2, with r2 (sufficiently) larger than n. Despite the wind be¬ 
ing nearly collisionless, Si is prevented from going through 
S2 by electromagnetic insta bilities, which gives rise to a CS. 
See lWaxman fe Loebl d200ll l for an estimate of their growth 
time. 

For an optically thin circumstellar medium, we con¬ 
firm nu merically the formation of a CS after SB (as 
found bv lKlein fc Chevalier! (I f 9 78| j and lEnsman fc Burrowsl 
lll992| yi with our ID-spherical radiation-hydrodynamics 
code, described in Section [I] In planar geometry, the r~ 2 
dilution factor is not prese nt, and no CS shou ld form, which 
agrees with the findings of lSapir et alj (|201lh . 


3 COLLISIONLESS SHOCKS BEFORE 

BREAKOUT FROM OPTICALLY THICK 

WINDS 

Let us now consider progenitors surrounded with optically 
thick winds. SB then starts in the wind at r ~ c/u s = 
/3s” 1 (IKatz et al.ll201(ih . 

We demonstrate below the central message of this pa¬ 
per: For some realistic optically thick circumstellar winds 
and for some shock velocities, the RDS does not survive long 
after leaving the hydrostatic core of the progenitor, and a 
CS forms at r > Ps ■ At r ^ /3s” 1 , radiation still plays a 
non-negligible role, and numerical simulations in the next 
Section show that such a shock resembles a ’decaying’ RDS 
containing a CS within its broad- width. 

Let us take r, ~ 10 13 cm (resp. ~ 10 11 cm) as or¬ 
ders of magnitude for radii of red supergiants (resp. Wolf- 
Rayets). The density profiles at r > r* in the optically 
thick winds are poorly known and may not be oc 1/r 2 . 
However, our results do not strongly depend on them. 
For the numerical simulations, we take p = M /47rw w r 2 
for a stellar mass loss rate M and wind velocity m w - 
Since r = k J r °° pdr, rbr ~ kM/? s / 47 tm w . Taking plausi¬ 
ble values for a dense wind of M « 5 • 10^ 4 Mgyr _1 and 
m w ~ lOkms - 1 (resp . 1000 km s^ 1 ) for RSG (resp. WR) 
progenitors dCrowtheri 120071 ). one can reach rhi/r, ~ 10 
for /3 S = 0.1, which is compatible with some interpreta¬ 
tions of SB observations dCampana et al . 2006; Wang et aid 


120071 : ISoderberg et al.1 120081 : Modiaz et al.l 2001 ). One may 
reach_£br/r*_^> 10 fo r sig nificantly larger M, see for exam¬ 
ple [Ofek^etj]!] d2010h and lSvirski et all d2012T ). Let us note 
that the following discussion is also valid for slower shocks 
with /3 S « 0.01. 

In this section, we simplify the problem to elucidate the 
essential physics. We assume that p = pi for r < r, (core). 
Since acceleration of wind shells depends only on the flux of 
radiation and not on the fluid density, we take for heuris¬ 
tic purposes p = p 2 for r* ^ r ^ rt, r (thick wind) with 
rbr r*, see Fig. Q] Al t = 0, the RDS leaves the core 
and enters the wind. Its width becomes ~ \ 2 c/u\ = A2 //3i, 
where Ai ,2 = l/«pi,2 is the photon mean free path and 
Mi the velocity reached by the shock after accelerating in 
the steep d e nsity g radient at the edge of t h e core -see 
e.g. ISaku rail digG dhlMatzner fc McKeaf dl999h . ISapir et ~al] 
d201lh . Ro fc Matznei ( 20131 4. In the thick wind, a CS ap- 
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Figure 1. Toy model of the problem, which captures the only 
relevant features for our analysis. The stellar core of radius 
r* = n (t = 0) and density pi = 1/reAi is surrounded with an op¬ 
tically thick wind with p2 < pi - The radiation-dominated shock, 
with velocity ~ u i and widths ~ \i/fh leaves the core at t = 0. 
Breakout occurs at r| )r 3> r». A collisionless shock forms when 
the shell at rj(t) catches up supersonically that at r 2 (t), which 
can happen at r -C r| )r , see text. 


pears if two shells catch up each other supersonically, i.e. 
the shell at lower radius rams into the other one with a rel¬ 
ative velocity (much) larger than the sound speed c s . The 
sh ock is collisionless rathe r than viscous, see calculations 
in IWaxman fc Loebl ihoQll h Micro-instabilities, such as the 
Weibel instability, mediate the formation of a CS. While this 
does not occur in a regular RDS, it can happen if fluid shells 
in the upstream cannot be accelerated to sufficiently large 
velocities by radiation to prevent this. 

Let us consider two shells <Si and 1S2 with initial radii 
ri(t = 0 ) = r* and r2(t = 0 ) = r* + d with d < A2// 3 i, see 
their initial locations in Fig. [T] With the above values, the 
mass swept up by the shock in the wind at r ^ rb r is negli¬ 
gible compared to the ejecta mass. Si then does not signifi¬ 
cantly slow down before SB and we can assume its velocity 
u 1 to remain constant. Let us denote by M2 the maximum 
velocity to which S2 can be accelerated by radiation only. In 
practice, c s -C ui because c s ~ 2 • 10 5 m s -1 yT^/lOOeV in 
the wind heated to T2. 5 i will then catch up 1S2 and create a 
CS before SB if (i) M2 < Mi — c s — mi and (ii) the radius by 
which they catch up is smaller than rb r : r* + < fb r . 

In the limiting case where no absorption of radiation oc¬ 
curs between Si and S2, all photons that have accelerated 
Si may accelerate S2. The velocity of a fluid shell is pro¬ 
portional to the integrated flux of radiation passing through 
it. Therefore, the maximum velocity reachable by S2 due to 
these photons does not exceed 

M 2 sph ^ Ml ( ' — ) . ( 1 ) 

\r, + d J 

The r ~ 2 dilution of radiation intensity is the main reason 
why Si may catch up S2. In practice, the actual velocity 
reached by S2 is larger or smaller than U2, sp h, depending 
on additional competing effects. First, photons pushing Si 
lose energy, and U2, sp h is likely to be overestimated in ®. 
Second, Ai increases in the expanding shocked core and ad¬ 
ditional radiation may leak out of it and accelerate both Si 
and S2. Third, the dense wind between the two shells may 
radiate through S2 part of its energy E em while being com¬ 
pressed, which further accelerates S2. The two first effects 
work in favour of Si catching up S2, while the third one has 
the opposite effect. Therefore, we estimate the latter one, so 


as to know if and when a collisionless shock can form. Since 
fj em ~ ff* +d Anr 2 ^u\ dr, an upper limit on the velocity of 
S 2 is 


M2 < Ml 


/ T * \ 2 K -Bern 

\r* + d/ c 47r(r* + d) 2 


This yields 


M2 ^ Ml 


1 + d 


i + Jl h+d 2 + ^- 
2A2 V 3 


( 2 ) 


where we have written lengths x in units of r* as x = x/r *. 
This estimate is likely to overestimate M 2 because the surface 
(oc r 2 (t) 2 ) of S 2 non-negligibly increases during its acceler¬ 
ation, and because not all kinetic energy will be radiated 
in reality. We perform numerical simulations in the next 
section for more accurate results. Radiation diffusing from 
the expanding core cannot prevent Si from catching up S 2 , 
because r 2 (t) Js n (t). However, a large extra acceleration 
might in principle make them meet beyond tv and prevent 
a CS from forming at r < rbr- We found this to have a 
negligible effect for relevant situations with our simulations. 

For an optically thick wind with density proportional to 
r -2 , d < A 2 //3i <C 1 is satisfied close to the hydrostatic core. 
Condition (i) then becomes Pi < 4 A 2 . Condition (ii) gives : 
Pi < 4A2[1 — l/2(rb r — 1)], which is not significantly more 
stringent than (i) for rb r ss 10. In practice, the main uncer¬ 
tainty lies in the factor ’4’, which is likely to be larger than 
this conservative estimate. From our numerical simulations, 
we find that approximately 


Pi < IOA 2 = 0.1 


f \ 

( r * 

{ M \ 

\ lOkm/s ) 

V 10 13 cm ) 

V5-10- 4 M o/yrj 


-1 


(3) 

for a r 2 wind (or equivalently, \ 2 IP 1 > r,/10), is a good 
overall estimate which does not noticeably depend either on 
the density profile of the progenitor, or on the sharpness of 
the transition between the core and the wind. Pi corresponds 
to the velocity of the shock when it enters the wind. For 
simplicity, we took in this analytical discussion a flat profile 
for the stellar envelope, but we verified numerically that 
our results hold for more realistic density profiles, such as 
p oc r -2 for a red supergiant envelope. 

Therefore, if the shock velocity does not exceed the 
value given by Eq. ©, the RDS does not survive the tran¬ 
sition from the core to the wind, because Si catches up S 2 
at r < rbr- A CS then forms in front of the expanding core. 
In other words, for a given wind density, and below a given 
shock velocity, the kinetic energy between Si and S 2 that 
can be radiated through S 2 is not sufficiently large to com¬ 
pensate for the dilution of photons due to the spherical ge¬ 
ometry of the problem. For progenitors with the above pa¬ 
rameters, this yields pi < 0.1, which is typical of RSG shock 
velocities. For such progenitors, a CS forms before SB. For 
some WRs, larger shock velocities can occur and a CS forms 
before SB only for larger values of X 2 C = r„): For example, 
the mildly relativistic shock of SN 2006aj does not satisfy 
Eq. © with M k, 3 • 10 -4 M@ yr _1 . 

More precisely, the scenario discussed here applies when 
the circumstellar wind is sufficiently dense to be optically 
thick, but not so thick that Inequality m is violated. For 
progenitors with steady winds (density p oc r~ 2 ), this 
typically corresponds to moderately thick winds with e.g. 
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rbr ~ 10 r*. Let us mention that in the hypothetical case 
of a progenitor surrounded by a thick wind with a density 
profile flatter than r~ 2 at r < rbr (e.g. due to variations in 
time of M before the explosion), this scenario can also be 
valid for significantly larger values of rb r /r*: It is valid as 
long as the condition /3i < IOA 2 is satisfied close to the core. 

Red supergiants or Wolf-Rayet stars with relatively high 
mass-loss rates prior to the explosion are good candidates : 
Moderately thick winds with rb r ~ 10 r* correspond to M 
in the range from a few x 10 -5 Mq yr -1 to ~ 10 -3 Mq yr -1 
for typical values of r«, u w and /3i for RSGs and WRs. For 
example, Type Ibc supernova SN 2008D/XRF 080109 may 
have been an event in which a CS is formed before SB : 
ISvirski fe Nakaii (120141 ^ suggest that SN 2008D is consistent 
with the explosion of a WR in a moderately thick wind, 
and that the progenitor underwent a steady but enhanced 
mass-loss during the l ast < 10 days. V e ry int erestingly, the 
parameters derived bv ISvirski fc Nakaii (l2014fl for SN 2008D 
(/3i ~ 0.25, u w « lOOOkms -1 , and Mm 2- lO _ 4 M 0 yr _1 
close to the star) marginally satisfy Inequality Q for r» « 
10 11 cm. This gives an additional and important reason to 
search for SN 2008D-like events in the future. The rate of 
such events depends on the likelyhood for a star to undergo 
enhanced mass-loss during the last few weeks or years pre¬ 
ceding the supernova. 

In contrast, progenitors of Type Iln supernovae are not 
expected to satisfy our scenario, since Eq. © , with typical 
values of M and u w for Type Iln, implies upper limits on /3i 
well below the actual shock velocities. For progenitors with 
significant mass-loss rates, the radiation-dominated shock 
should always survive the transition from the core to the op¬ 
tically thick wind, and the picture would then be a conven¬ 
tional one : In such cases, photons are supplied by the imme¬ 
diate downstream of the RDS in the thick part of the wind. 
For example, Type Iln supernova SN 2010)1 does not satisfy 
Eq- © because of the large mass-lo ss rate M ~ 1 Mq yr 1 
dZhang et al.ll2012l ; lOfek et al.ll2014l l. 

In stellar cores, RDS are stable because A is sufficiently 
small to prevent conditions similar to those of Eq. © from 
being met. 


4 NUMERICAL SIMULATIONS 

We confirm numerically the predictions of Section [3] with 
our Eulerian ID-spherical radiation-hydrodynamics code. 
The fluid is assumed to be fully ionized. The code is two- 
temperature, i.e. we assume proton and electron tempera¬ 
tures to be equal. The results presented below are not af¬ 
fected by this assumption. We use a gray frequency aver¬ 
age for the radiation, and represent it by its internal energy 
E ra d with characteristic temperature T ra d = (cE ra d/40-) 1 / 4 , 
where a is the Stefan-Boltzmann constant. At each time 
step, the radiation transport is solved using a ’square- 
root’ (M orell 12000 1 flux-limited diffusion approximation, 
and the opacity k is assumed to be dominated by Thom¬ 
son scattering. For the transfer of energy between fluid 
and radiation, we take into account Co mpton cooling and 
brems strahlung, using the formulas of IChevalier fe Kleinl 
d 19791) . Underl ying assumptions for a n equivalent code can 
be found in IChevalier fe Kleinl dl979l ). 

Our main result does not strongly depend on the density 


and temperature profile of the hydrostati c core, and we take 
as an example those for the RSG used in IChevalier fe Kleinl 
dl979i l. We choose initial conditions such that the shock ve¬ 
locity reaches /? s ft; 0.1. We use winds with density pro¬ 
files oc r ~ 2 . We take here M = 5 ■ 10~ 4 MQyr _1 and 
u w = 10kms _1 , so rb r ft; 3. A CS appears before SB, at 
f < 1.5. At the CS formation time, photons still have not 
started to escape from the optically thick material, con¬ 
trary to e.g. expectations for breakout from a stellar sur¬ 
face. Fig. [2] shows the CS at f ft; 1.6, near the downstream 
of the radiation-dominated transition (remains of the ini¬ 
tial RDS). See caption of Fig. [2] for details. It appears as a 
growing discontinuity in the smoother velocity profile. The 
radiation-dominated transition extends to radii larger than 
shown in Fig. [2] At such an early time, the CS downstream 
temperature is only ~ 1 keV because the radiation still pro¬ 
vides most of the fluid acceleration in its upstream, but we 
find the discontinuity in u to grow and the CS processes 
a significantly larger fraction of pu 2 at larger r. For other 
parameter values (smaller /3 S for M, u w fixed), radiation 
plays a smaller role, allowing the CS to emit in hard X- 
rays (> 10s keV) before its photons break out. The formula 
/3 S < IOA 2 and conclusions of the previous section have been 
verified by scanning the parameter ranges. We find that for 
significantly larger M = (1, 5) ■ 10~ 3 Mq yr -1 and the same 
core profile, the RDS survives the transition to the wind. 
Details of the hydrostatic core are not found to be very im¬ 
portant, and we test the WR case by rescaling, as a first 
approximation, the above profile to r* = 10 11 cm and larger 
densities. We vary /3 S by slightly changing the explosion en¬ 
ergy. For M — 5- 10 _4 MQyr^ 1 and w w = 1000kms -1 , we 
find that for /3 S > 0.15, no CS appears before SB, whereas 
they do appear before SB for /3 S < 0.1. Varying M and u w , 
also corroborates the picture announced in the previous sec¬ 
tion. 

We verified numerically that once the CS is formed 
(even when this happens at r < rbr), it survives to r > rbr 
in winds with p oc r -2 : If the shock does not sufficiently 
slow down at r $ rb r and if the wind density does not in¬ 
crease with r, no process reduces the difference in velocities 
between the CS immediate downstream and immediate up¬ 
stream down to a subsonic value. In the simulations, this 
difference in velocities is, on the contrary, found to grow at 
r < r br . 


5 PARTICLE ACCELERATION 

Assuming conservatively a magnetic field strength at 
the C S similar to t h at at the stellar surface, B s ~ 
10 G llBarvainis et alj Il987i l. and wind densities p ~ 
IQ-ii (-9) gem -3 , the CS is super-Alfvenic. Once it is 
formed, CR acceleration may start. Coulomb losses for 
suprathermal particles are sufficiently small here and do 
not prevent them from entering diffusive shock accelera¬ 
tion and being accelerated. However, for WRs with M > 
10 _ 3 MQyr _1 , such losses start to inhibit CR accelera - 
tion before SB. Some findings of W axman fe Loei] (1200 il l 
and iKatz et alj d201ll) can be transposed to our study, 
yet we deal here with a shock propagating in denser re- 
gion s of the wind. Assuming Bohm diffusion for CRs at the 
CS dReville fc Belli 120131 : ICaprioli fe Spitkovskvll 20 lJ ). and 
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Figure 2. Simulation of a RSG explosion in an optically thick wind before shock breakout, see text for parameters. Zoom around the 
downstream of the radiation-dominated transition of the shock, where a discontinuity (CS) can be seen around r £3 5.28 • 10 11 m, at 
1.6 r* and « r*b r /2. Left panel: Most of the fluid acceleration d(pu)/dt is due to the radiation, except in a thin zone at the CS, where 
the fluid contribution — Vp dominates. It coincides with the sharp discontinuity in u. Right panel: A dense shell appears in p, in the 
downstream of the CS. Peaks in p and in the electron temperature T e appear in the CS immediate downstream. T ra d follows T e except 
in the peak region, where it remains flat —Frad 1S smooth. 


equal dwell times in the downstream and the upstream, one 
finds a typical acceleration time 

r CR ~ 7T~^ o ~ 30 S ( FFFfr ) ( ) 7^ ' (4) 


3 eB s ui 


10 G 


for protons. This time can be optimistic when the disconti¬ 
nuity in velocity at the shock is still small due to smooth¬ 
ing by radiation. However, magnetic field amplification at 
the shock due to the non-resonant hybrid (NRH) instabil¬ 
ity (iBelll 120041 1 plays a role in the opposite direction by 
di minishing tqr and ther eby facilitating CR acceleration, 
see iGiacinti fe Belli (l2014Tl for a detailed study. Magnetic 
field amplification is (constantly) driven by the escape of 
the highes t energy CRs in t he upstream of the collisionless 
shock, see iBell et ahl (120131 1. For the ranges of parameter 
values that are relevant here, the typical growth time of the 
NRH instability is smaller than the damping time of the 
turbulence by the radiation field, which energy density is 
U rad ~ puf Therefore, magnetic field amplification should 
occur in such conditions. In the upstream of the CS, a tur¬ 
bulent fluid parcel with velocity Ut suffers momentum losses 
due to radiation (second order Fermi for photons). From the 
momentum equation of the fluid parcel, one can deduce the 
typical damping time of the turbulence : 


7"damp — 


U t 


du t /dt 


KpUsUt 


> 


KjpUi 


(5) 


The size of the discontinuity in velocity at the CS may be 
written as Am = |/u s , where 0 < / ^ 1 and / = 1 is the 
limiting case where no radiation accelerates the upstream 
of the CS. The growth rate of the fastest growing mode 
of the NRH instability is equal to 7 ma x = 0.5j'cr \JPo/p, 
where j'cr — 0. 03 pf 2 vfe/EcR i s the CR current density 
which drives it (iBell et al.l l2013l l. The instability growth 
time, tnrh ~ 57^1x1 is then small compared to Tdamp : 

tnrh < 10 ^ ( F C r \ / _ J _\ ~ 2 

~ 0.03 c 2 e % /po/ 2 ~ ' V 10TeV / V°- 05 / 


_ M_ _\ 1/2 / 

5 ■ 10~ 4 M Q /yr ) Vl° km /s 

numerically for a wind with p oc r~' 

Let us note that lOTeV energies are reachable before 
breakout because tor <C (fb r — r*)/u s ~ several hours (resp. 
minutes) for RSG (resp. WR) progenitors with /3 S = 0.1 and 
rb r ~ 10. For such RSGs, rcR(at lOTeV) is smaller than 
energy loss times from pion production through inelastic pp 
and py collisions. The typical life time of a CR proton due 
to pp collisions, r pp ~ m p /0.2cp<j pp , is 



-1/2 


(io^r (6 > 


: 4 min 


f U-w \ 

( r y I 

f M \ 

\ 10km/s ) 

Vl0 13 cm/ ' 

V 5 • 10~ 4 M 0 /yr J 


(7) 


The background ~ 10 eV photons in the thick wind are 
not sufficiently energetic to trigger pion production through 
inelastic P 7 scattering. For lOTeV CRs, > lOkeV photons 
are required to exceed the threshold for pion production. 
Photons with such energies can be produced by the radiative 
CS. However, the number de nsity of target ph otons n 7 must 
be much less than pu^/hv dKatz et al.l 1201 il l. We find for 
the typical life time of a CR proton due to pq collisions, 

T p7 CS I / 0.2 cM-, ct p7 . 


r P 7 ^ 2 min 


( \ 

( r V l 

( M \ 

\10km/s J 

Vl0 13 cm7 ' 

\ 5 • 10 -4 M 0 /yr J 


-1 
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(JLY 2 ( Ec R 

VO.l/ \ 10 TeV 



( 8 ) 


e* pair creation due to py interactions does not yield a 
stronger constraint. 

In the case of Wolf-Rayet progenitors with the above 
parameters, t PP:P 7 > 3 s. Consequently, TeV energies may 
be reached for WRs. 
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6 OBSERVATIONAL CONSEQUENCES 

A discussion on progenitors for which the CS is expected 
to form before SB may be found at the end of Section [3J 
from Eq. ©• We now describe the two main observational 
consequences : X-ray flashes and high-energy neutrinos. 

For fast shocks (/ 3 S > 0.1 — 0.2), RDS start to depart 
from thermal equ ilibrium, which may produce X RFs in asso¬ 
ciation with SBs (IWeaveijll976i : [Sapir et al.ll201ll f. For slower 
shocks, UV photons are typically expected. We predict that 
even for lower /3 S , photons with energies > (1 — 10s) keV can 
be emitted from the beginning of SB, but for a different rea¬ 
son: This happens when the radiative CS forms significantly 
before SB, such as for RSGs (or WRs with ’slow’ shocks) sur¬ 
rounded with optically thick winds. It would heat the plasma 
at r <C rbr to temperatures higher than expected for a ’slow’ 
RDS (~ 10 —100 eV) and load the thick wind with (hard) X- 
rays. This would result in a flash at breakout that both con¬ 
tains hard photons reflecting the presence of the hot down¬ 
stream of the CS, and softer photons (notably LTV) from the 
remains of the former RDS. The energy radiated at breakout 
is typically ~ 10 ^^_erg depen ding on the tested pro geni- 
tors -see also iKatz et all d2012l l an d lSapir et all (120131 ). We 
find that, from the beginning of SB, the fraction of the en¬ 
ergy emitted in X-rays is already roughly comparable with 
that in softer photons: From ss 10% to more than a half, 
with the largest fractions also corresponding to the high¬ 
est maximum X-ray energies. The X-ray flux rises abruptly 
on a time scale ~ rv/c « 30 -), and then decays 

more slowly due to the persistence of inverse Compton on 
background photons in a wind with p oc r -2 . Observations 
in different energy bands will be needed to distinguish be¬ 
tween progenitors following the scenario presented here, and 
those for which the CS only starts to form during breakout, 
as suggested bv lKatz et al.1 d201l] ~). The production history 
of CRs at t T§> /3a” 1 should also leave imprints in the spec¬ 
trum at SB : Secondary 7 -rays, notably from 7 r° decay, are 
injected in the wind and partly reprocessed to lower ener¬ 
gies through e ± cascades on the large photon background 
(7 + 7b -» e + + e~). 

The fact that outer layers of the thick wind at r < rbr 
may have been mostly accelerated by the CS implies that 
the energy radiated at SB may be 47 rr 2 r c 2 /3 S / k. This may 
ease the tension between radiated ene rgy and duration of 
the emission for XRF 080109 (see e.g. ISapir et al.l (l2013h l . 
provided this event corresponds to SB from such an optically 
thick wind. The relatively low photon flux at breakout would 
be consistent with more energy being in the thermal plasma 
behind a CS, as expected in our scenario. 

This work also provides a new technique to access infor¬ 
mation on SN progenitors inside thick winds, such as the ra¬ 
dius of the stellar core r», and the density profile at r > flif 1 . 
By detecting secondary > 100 GeV — 1 TeV neutrinos (from 
notably n ± decay) before the first photons from breakout, 
one will improve our knowledge of the still poorly under¬ 
stood late stages of massive star evolution. The time inter¬ 
val between the arrival of the first neutrinos and photons 
is Af „ 7 fts (rbr — ?'»)(Ar 1 — l)/c ~ 8 hr (resp. 5 min) for 
RSGs (resp. WRs) with the above parameters, rb r = 10 
and /3 S = 0.1. Assuming that 5 % of the energy processed 
by the shock is channelled into CRs, we typically find for 
a source at distance l, and a processed mass between r* 


and rbr of « 10 ~ 5 M@, that ~ 10 3 (3kpc/I ) 2 neutrinos with 
~ TeV energies would be detectable before SB by IceCube or 
KM3NeT. One could record a few of such neutrinos for an 
event in the Magellanic Clouds. The low rate of such super¬ 
novae within ~ 100 kpc from Earth is the main limitation to 
the detection of these neutrinos with a km 3 observatory. For 
example, the rate of 2008D-like supernovae in our Galaxy 
should be at most 1/1000 yr. 

A supernova detected in neutrinos in the pre-shock 
breakout phase will generally yield more neutrinos in the 
post-shock breakout phase, except if the density of the pro¬ 
genitor wind suddenly falls sharply just beyond rbr- These 
latter neutrinos will be detected after SB has started. For ex¬ 
ample, in a steady wind with density p oc r -2 , the mass pro¬ 
cessed in the post-shock breakout phase by a shock travelling 
between ry, r and r > rb r , is ~ M(r — ry r )/uw This amount 
quickly exceeds that processed i n the pre- s hock breakout 
phase (~ M(ry r — r*)/u w ). iMurase et al.l d201ll l studied 
in detail the post-shock breakout emission of neutrinos for 
shocks interacting with dense circumstellar material, such 
as shells. 


7 CONCLUSIONS 


During a core-collapse supernova, a radiation-dominated 
shock propagates through the progenitor star. If the sur¬ 
rounding wind is optically thin, this shock stalls when it 
reaches the outer layers of the stellar core. In the upstream, 
the circumstellar material is then accelerated by escaping 
photons from shock breakout to a velocity roughly oc r~ 2 , 
where r is the distance to the centre of the progenitor. The 
outer layers of the shocked core ram supersonically into these 
slower layers of the wind, and a collisionless shock is ex¬ 
pected to form during or on th e time scale of supernov a 
shock breakout. See, for ex ample, Chevalier fe Kleinf d l979h . 
lEnsman fc Burrows! (ll992T ) and IWaxman fe LoebT i 200lh for 
detailed studies. 

In the present paper, we have investigated the case of 
supernovae occuring in thick winds. In this case, the forma¬ 
tion of a CS should also occur no later than during or on 
the time scale of shock breakout —from the ’outer’ layers of 
the optically thick part of the wind. 

We have demonstrated here that for some 
astrophysically-relevant progenitors surrounded with 
thick winds, a collisionless shock forms well before break¬ 
out, providing new ways to study invisible layers of their 
winds and to constrain stellar evolution theories. In such 
cases, the RDS has been found to stall when entering the 
optically thick part of the wind, notably because of shock 
curvature. Photons are then mostly supplied by the shock 
propagating in the core, and the wind is not sufficiently 
dense to compensate for the r -2 dilution of photons in 
the wind. On the contrary, for progenitors where the RDS 
survives the transition from the core to the wind, such as 
for type Iln supernovae, photons are mostly supplied by the 
immediate downstream of the shock in the wind. 

We have discussed, in Section[5] the onset of particle ac¬ 
celeration at the CS. For example, we predict that for some 
red supergiants surrounded with thick winds, a fraction of 
secondary high-energy neutrinos from CRs can arrive ~ 10 
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hours before photons from shock breakout, and more neu¬ 
trinos are produced later in the post-shock breakout phase. 

We find that the CS forms after the RDS exits the 
core, and before breakout, for progenitors with shock ve¬ 
locities < 0 . 1 c ( T o^)( 5 , i 0 -4 M Me/yr )~ 1 ( T onTiir). where “w, 
M and r* respectively denote the wind velocity, mass- 
loss rate and radius of the hydrostatic core. The wind has 
to be sufficiently dense to be optically thick but not ex¬ 
cessively. For progenitors with steady winds (p oc r _2 L 
this corresponds to moderately thick winds, where e.gjj 
rbr ~ 10 r*. Progenitors of Type Iln supernovae are ex¬ 
pected to have too dense winds to form CS when the RDS 
leave their cores. However, Wolf-Rayet stars or red super¬ 
giants with either dense winds or enhanced mass-loss prior to 
the explosion are better candidates. For example, Type Ibc 
su pernova SN 2008D/XR F 080109 has been interpreted 
by lSvirski fe Nakad l|2014l l as the explosion in a moderately 
thick wind of a WR star, undergoing an enhanced mass- 
loss during its last < 10 days . Inter estingly, the parameters 
inferred by ISvirski fe Nakad (120141 ') for SN 2008D make it 
marginally consistent with the above condition. This is an¬ 
other important motivation to search for similar events. In 
the future, one can notably use them to study the forma¬ 
tion times of collisionless shocks with respect to the photon 
flashes at breakout. 

More generally, supernovae occurring in dense winds are 
promising targets for multi-messenger studies. The detection 
of their UVs, X-rays, 7 -rays and TeV neutrinos will allow 
one to test a wide variety of physical and astrophysical phe¬ 
nomena in extreme conditions, such as particle acceleration, 
magnetic field amplification and shock physics. 

Finally, studying CR acceleration in dense winds is im¬ 
portant, because it should lead to a bet ter understanding o f 
the knee in the C R spectrum, see e.g . Svcshnik oval (I2003T) . 
iBell et ahl (120131 ') , iMurase et all (1201 4l ) . 
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